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ABSTRACT 



m ■ 
o , 

■ We employ quantum mechanical principles in the computability exploration of the class of classically noncom- 
putable Hilbert's tenth problem which is equivalent to the Turing halting problem in Computer Science. The 
Quantum Adiabatic Theorem enables us to establish a connection between the solution for this class of problems 
and the asymptotic behaviour of solutions of a particular type of time-dependent Schrodinger equations. We 
then present some preliminary numerical simulation results for the quantum adiabatic processes corresponding 
to various Diophantine equations. 



1. HILBERT'S TENTH PROBLEM 

One of the well-known and fascinating mathematical problems is about the existence of integer solutions of 
polynomial equations of the type 



g: {x + lf + {y + lf -{z + lf = 0. 

According to Fermat's last theorem, which has only been proved very recently, the above equation does not have 
any integer solution for x, y, and z. However, although it took hundreds of years for mathematicians to finally 
come up with such a proof, this proof may not be valid for similar equations such as 



{x + \f + {y + \f -{z + \f + xyz = 0. 



^ : 

^ ■ Many important mathematical problems, such as Goldbach's conjecture or the distribution of zeroes of Riemann 
Zeta function, can be cast into equivalent problems of whether related polynomial equations with integer coeffi- 
■ cients have integer solutions or not. Because of this, David Hilbert in his famous list of mathematical problems 
. ^ in 1900 included as the tenth problem^ the challenge: 

><: . 

, Given any polynomial equation with any number of unknowns and with integer coefficients: To 

■ - - ' devise a universal process according to which it can be determined by a finite number of operations 

whether the equation has integer solutions. 



Such equations are known as Diophantine equations. Never was it anticipated that this tenth problem is ul- 
timately equivalent (as shown by Davis, Putnam, Robinson and Matiyasevich, see^) to the halting problem of 
Turing machines of more than 30 years later. On this equivalence basis it has been concluded that Hilbert's 
tenth is not computable: there is no single universal process to determine the existence of integer solution or lack 
of it for arbitrarily given Diophantine equations in as far as there is no single universal machine to determine 
the halting or not of arbitrarily given Turing machine (which starts with some arbitrarily given input). 

Thus, we would have to consider anew each different case of Diophantine equation. 

In spite of this widely accepted result, we have proposed a quantum algorithm? for solving Hilbert's tenth 
problem. In the next section we briefly summarise the algorithm. We then next present some preliminary 
numerical results from the simulations of quantum processes for some very simple Diophantine equations as a 
concrete "proof of concept" . 
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2. A QUANTUM ALGORITHM 

2.1. General oracle for Hilbert's tenth problem 

It suffices to consider only nonnegative solutions of a Diophantine equation. Let us consider the example 

[x + if + {y + if ~{z + If + cxyz = 0, c&Z, 
with unknowns x, y, and z. Starting from the observation that if we can construct the hamiltonian 

Hp = ((4 ax + 1)^ + [alay + if - {ala, + if + c{aia^){alay){ala^))^ , 

which has a spectrum bounded from below in fact, and if we can obtain the corresponding ground state (of least 
energy) then we can solve the tenth problem! 

The ground state \g) of the hamiltonian so constructed has the properties, for some {nx,ny,nz), 
Nj\9) = nj\g), 

Hp\g) = {{n:o + lf + {ny + lf -{n^ + lf + cna:nynzf\g)=Eg\g). 

Thus a projective measurement of the energy Eg of the ground state \g) will yield the answer for the decision 
problem: The corresponding Diophantine equation has at least one integer solution if and only if Eg = 0, and 
has not otherwise. (If c = in our example, we know that Eg > from the Fermat's last theorem.) 

Prom the above observation, our ground-state oracle is thus clear: 

1. Given a Diophantine equation with K unknowns x's 

D{xi,---,Xk) = 0, (1) 
we need to simulate on some appropriate Fock space the quantum hamiltonian 

Hp = (^D{a\au---,a^K(iK)y ■ (2) 

2. Measurement results of appropriate observables in the ground state will provide the answer for our decision 
problem. 

One way, which is by no nu^an the only way. to obtain the ground state is guaranteed by the quantum 
adiabatic theorem"' which we will exploit in the next section. 

2.2. Quantum Adiabatic Algorithm 

In the adiabatic approach,'"' one starts with a hamiltonian Hj whose ground state \gi) is readily achievable. 
Then one forms a "slowly" varying hamiltonian which interpolates between Hj and Hp in the time interval 

te[o,T] 

n{t/T) = {l-t/T)Hi + {t/T)Hp, 

= Hj + {t/T)W. (3) 

We will adopt this approach with the proposed (universal) initial Hamiltonian 

K 

Hi = J2(4-^:)i(^i-^i) (4) 

which admits as the ground state the coherence state 

\gi) = \ai) igi ...\aK), (5) 



where = e 2!"^' ^^^^ 

Provided the conditions of the adiabatic theorem'' are observed, the initial ground state will evolve into our 
desirable ground state \g) up to a phase: 

^lirn^ Texp I -iT^'H(r)dr I = e'^\g), (6) 

where T is the time-ordering operator. 

So our problem now is to solve the time-dependent Schrodinger equation for t G [0, T] 

idtm = n{t/T)i^{t), [0,T], (7) 

with the initial state being the coherence state (5). 

We can analytically show in general the two crucial results below^: 

• The ground state of H{t/T) is non-degenerate for t G (0,r). As the minimum energy gap between the 
ground state and the first excited state is non-zero, it takes only a finite time T for the adiabatic process, 
as asserted by the quantum adiabatic theorem, to generate a state which has a high probability of being 
the ground state of Hp. 



The probability of the state at time T in some number state, |(V'(T')|no)| , is greater than 1/2 iff \no) is 
the ground state of Hp. That is, we only need to solve the equation for increasing T until this majority 
condition is satisfied in order to identify the ground state. 

The proofs of these important results will be available elsewhere. 

3. SIMULATION TECHNICALITIES 

• We solve the Schrodinger equation numerically in some finitely truncated Fock space large enough to 
approximate the initial coherence state to an arbitrarily given accuracy. That is, a truncation m is chosen 
such that 



a,;m) = e-^l"'l V^|n) (8) 



n=0 



has a norm less than one by some chosen e, 1 1 — (a^; m\ai] m) \ < e. 

At each time step St, and up to 0{5t^), 

i){t + St) = exp{-in{t/T)St}ip{t) = (1 - in{t/T)St - l/2H^{t/T)St'^)i}{t), (9) 

thus there are maximally only two creation operators, aja], and we can explore this fact to explore the 
infinite Fock space by increasing the size of the truncated Fock space by two at every time step. 

We employ the unitary solver, with conjugate gradient method, 

' 1 + ^H{t/T)St^^ ' ^ ' 

which approximates exp{—iH{t/T)St} to second order in St. This solver, being unitary, preserves the norm 
of the state vector in the evolution of time even though the size of the underlying truncated Fock space is 
allowed to increase with time if necessary. 

The time step St at time f is a function of time and is chosen in such a way that halving this St{t) only 
results in 0{St^) correction. 
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Figure 1. Corresponding equation xy + x + 4y ~ 11 = 0; with intermediate evolution time T. 



4. SIMULATION PARAMETERS AND RESULTS 

4.1. Equation xy + x + 4y — 11 = 

This is an equation resulted from the factoring of the number 15 into two prime factors. Because of the nature 
of the problem we can fix the truncated size of our Fock space to be {m^, my) = (9, 9) in order to have the norm 
of the state vector less than unity by an amount e = 10~^ at all times. In this and in all other simulations below 
we choose ax — cny — (2.0, 0.0). 

In Fig. 1 we plot the magnitude square of maximum component (in terms of the Fock states) of the state 
vector as a function of the total evolution time T in some arbitrary unit. Closed to below T ^ 1500, the 
maximum probability components of \^(T)) are dominant by some states, denoted by (blue) star and (green) 
triangle symbols, none of which are the ground state. In fact, they are, respectively, the first degenerate excited 
states \nx,ny) = |4, 1) and |3, 1) ol Hp. 

With the increase in T in Fig. 2, one of the Fock state, |1, 2) denoted by (red) box symbol, has probability 
greater than ^ which is our criterion for being identified as the ground state. We mark this regime as the quantum 
adiabatic regime when the ground state wins the battle for dominance. From the ground state so identified we 
can infer that our Diophantine equation has one solution in this domain. 

4.2. Equation x + 20 = 

We consider this extremely simple equation as an example which has no solution in the positive integers. The 
simulation parameters are as in the previous, except that the truncated size of our Fock space, starting with size 
8, is now allowed to vary with time in order to simulate the exploration of the whole infinite space. 

In Fig. 3 we plot the probabilities of the dominant components as a function of T. Below T ^ 50 none of the 
two components is greater than one-half, and in fact the first excited state, |1) denoted by (blue) triangle symbol, 
clearly dominates in this regime. Eventually we enter the quantum adiabatic regime upon when the (red) box 
symbol rises over the one-half mark; indeed it corresponds to the Fock state |0) which is the true ground state 
and which implies that our original Diophantine equation has no integer solution at all. 

4.3. Equation x - 20 = 

We now consider a simple example which nevertheless has all the interesting ingredients typical for a general 
simulation of our quantum algorithm. Here we choose e = 10^"^ and the initial Fock space has only up to 
|n) = |14), which does not include the true ground state of Hp. This is typical in our simulations since we in 
general would not be able to tell in advance whether our initial Fock spaces do contain the true ground states or 
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Figure 3. Corresponding equation x + 20 = 0. 
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Figure 4. Corresponding equation x — 20 — 0. 
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Figure 5. The expectation value of occupation number, (A'^), as a function of T. 



not. Generally, they do not. Even so, our strategy of allowing the expansion in the size of truncated Fock space 
in time has enabled the true ground state to be found and identified. 

In this example, the state |20), which is not included in the initial truncated Fock space, is eventually reached 
and identified as the ground state as shown as red boxes in Fig. 4. Blue triangles and green stars are corresponding 
to the first two excited states |I9) and |2I), which are degenerate eigenstates of Hp. 

Note that these competing pretenders somehow have unexpectedly probabilities greater than one-half (around 
T ~ 90), contrary to our analytical result that only the ground state can have probability rising above one- half! 
We think that this is only some artefact of finite-size time steps St, and expect that it would go away once we 
employ a more sophisticated method for solving the Schrodinger equation. Work is in progress to systematically 
extrapolate to zero-size time steps to confirm the removal of this type of finite-size effects. 

Figs. 5 and 6 depict the expectation values of occupation number and of energy as functions of T. They are, 
respectively, approaching 20 and zero (which signify the fact that our equation has an integer solution, namely 
20). 
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Figure 6. The expectation value of energy, (Hp), as a function of T. 



5. CONCLUDING REMARKS 

We here consider the issue of computabiUty in principle, not that of computational complexity. The simulation 
results support the idea of implementation and execution of our quantum algorithm with a physical process^: 

• Run the physical process (corresponding to the Diophantine equation in consideration) for some time T. 

• Repeat the process at this T to obtain the statistics through measurements. 

• If none of the measurement outcomes exhibits probability of more than 1/2, ram up the T and go back to 
the step above. 

• Eventually, at some sufficiently large T, the measurement state with more-than-even probability can thus 
be identified as the ground state, terminating our physical implementation and execution. 

• Substituting the quantum numbers for the now identified ground state enables us to see ii Eg — and thus 
whether the Diophantine equation has a solution. 
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